QCD corrections to J/ip polarization of hadronproduction at Tevatron and LHC 
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The next to leading order (NLO) QCD corrections to J/ip polarization of hadronproduction 
at Tevatron and LHC are calculated. The results show that the J/ip polarization is extremely 
changed from more transversal polarization at leading order (LO) into more longitudinal polarization 
at NLO. Although it gives more longitudinal polarization than the recent experimental result on 
the J/ip polarization at Tevatron. It sheds light on the solution to the large discrepancy of J/ip 
polarization between theoretical predication and experimental measurement, and suggests that the 
next important step is to calculate the NLO correction for color octet state J/tp^ hadronproduction. 
Our calculations are performed in two ways where the polarizations are summed analytically or not, 
and they are checked with each other. It also gives a K factor for total cross section (ratio of NLO 
to LO) of about 2 and shows that the NLO corrections boost the J/ip production for about 2 order 
of magnitude in high transverse momentum p t region of J/ip, which confirms the calculation by 
Campbell, Maltoni and Tramontano. 
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The study of J/tp production on various experiments is 
a very interesting topic since its discovery in 1974. It is a 
good place to probe both perturbative and nonperturba- 
tive aspects of QCD dynamics. To describe the huge dis- 
crepancy of the high-p t J/ip production between the the- 
oretical calculation and the experimental measurement, 
color-octet mechanism [l| was proposed based on the non- 
relativistic QCD(NRQCD)[2|. The factorization formal- 
ism of NRQCD provides a theoretical framework to the 
treatment of heavy-quarkonium production. It allows 
consistent theoretical prediction to be made and to be 
improved perturbatively in the QCD coupling constant 
a s and the heavy-quark relative velocity v. Although it 
seems to show qualitative agreements with experimental 
data, there are certain difficulties in the quantitative es- 
timate in NRQCD for J/ip and ip' photoproduction at 
the DESY ep collider HERA @,|||, J/ip(ip') polarization 
of hadronproduction at the Fermilab Tevatron, and J / ip 
production in B-factories. 

There are a few examples shown that NLO corrections 
are quite large and it is difficult to obtain agreement be- 
tween the experimental results and leading order the- 
oretical predictions for J/ip production. It was found 
that the current experimental results on inelastic J/ip 
photoproduction are adequately described by the color 
singlet channel alone once higher-order QCD corrections 
are included [|. Ref. @ found that the DELPHI @ data 
evidently favor the NRQCD formalism for J/ip produc- 
tion 77 — > J/ipX, but rather the color-singlet mecha- 
nism. And it was also found in ref. that the QCD 
higher order process 77 — ► J/ipcc gives the same order 
and even larger contribution at high p t than the leading 
order color singlet processes. In ref. [g], at NLO the pro- 
cess eg — > J/ipc where the initial c quark is the intrinsic c 
quark from proton at Tevatron, gives larger contribution 
at high p t than the leading order color singlet processes. 
The large discrepancies found in the single and double 



charmonium production in e + e~ annihilation at B fac- 
tories between LO theoretical predictions 0, [l(| and ex- 
perimental results [Til . Il2| were studied in many work. 
It seems that it may be resolved by including higher or- 
der correction: NLO QCD and relativistic corrections 

anus. 

Based on NRQCD, the LO calculation predicts a siz- 
able transverse polarization for J/ip at high pt[15j while 
the measurement at Fermilab Tevatron [16[ gives almost 
unpolarized result. Recently, NLO QCD corrections to 
J/ip hadronproduction have been calculated by Camp- 
bell, Maltoni and Tramontano 17] . The results show that 
the total cross section is boosted by a factor of about 
2 and the J/ip transverse momentum p t distribution is 
enhanced more and more as pt becomes larger. A real 
correction process gg — > J/ipcc at NLO was calculated 
by Artoisenet, Lansberg and Maltoni [lf|. It gives sizable 
contribution to pt distribution of J/ip at high p t region, 
and it alone gives almost unpolarized result. A s-channcl 
treatment to J/ip hadronproduction gives longitudinal 
polarization by H. Haberzettl and J. P. Lansberg [l9j 
Therefore it is very interesting to know the result of J / ip 
polarization when NLO QCD corrections are included. 
In this letter, we calculate the NLO QCD corrections to 
the J/ip polarization in hadronproduction at Tevatron 
and LHC. In the calculation, we use our Feynman Dia- 
gram Calculation package (FDC)[20] with newly added 
part of a complete set of method to calculate tensor and 
scalar integrals in dimensional regularization, which was 
used in our previous work[l4j|. 

For LO process g{pi) + g(p 2 ) -> J/ip(pz) + g{pi), by 
using the NRQCD factorization formalism, the partonic 
cross section is expressed as 



da u 
dt 
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with 
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where i? s (0) is the radial wave function at the origin of 
J/ip and the approximation Mj/^ = 2m c is taken. The 
LO total cross section is obtained by convoluting the par- 
tonic cross section with the parton distribution function 
(PDF) G g (x,fj,f) in the proton: 



dx 1 dx 2 G g (x 1 , [if ) G g (x 2 , Hf ) a 



~B 



(2) 



where /i/ is the factorization scale. In the following, a 
represents corresponding partonic cross section. 

The NLO contribution to the process can be separated 
into the virtual corrections, arising from loop diagrams, 
and the real corrections, arising from the radiation of a 
real gluon or a light (anti)quark or the final-state gluon 
splits into light quark-antiquark pairs. There are UV, 
IR and Coulomb singularities in the calculation of the 
virtual corrections. UV-divergences from self-energy and 
triangle diagrams are canceled upon the renormalization 
of the QCD gauge coupling constant, the charm quark 
mass and field, and the gluon field. Here we adopt same 
renormalization scheme as ref. [21). The renormaliza- 
tion constant of charm quark mass Z m and field Z 2 , 
and gluon field Z3 are defined in the on-mass-shell(OS) 
scheme while that of QCD gauge coupling Z g is defined 
in the modified-minimal-subtraction(MS) scheme: 
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Where je is Euler's constant, (3q — ^-Ca — ^Tprif is 
the one-loop coefficient of the QCD beta function and 
n/ is the number of active quark flavors. There are three 
massless light quarks u,d,s, so nf—3. In SU(3) C , color 
factors are given by Tp = \,Cp = |, Ca — 3. And (i r is 
the renormalization scale. 

After having fixed the renormalization scheme, there 
are 129 NLO diagrams , including counter-term diagrams, 
which are shown in Fig. [TJ divided into 8 groups. Dia- 
grams of group (e) that has a virtual gluon line connected 
with the quark pair lead to Coulomb singularity ~ ir 2 /v, 
which can be isolated by introducing a small relative ve- 
locity v = \p c — Pc\ and mapped into the cc wave function. 
By adding all diagrams together, the virtual corrections 
to the differential cross section can be connected to the 
virtual amplitude as 



d^ 
dt 



cx 2Re(M B M v *) 



(4) 




FIG. 1: One- loop diagrams for gg — > J/ipg. Group (a) and 
(b) are counter-term diagrams of the quark-gluon vertex and 
corresponding loop diagrams, Group (c) are the quark self- 
energy diagrams and corresponding counter-term ones. More 
diagrams can be obtained by permutation of external gluons. 



where M B is the amplitude at LO, and M v is the renor- 
malized amplitude at NLO and is UV and Coulomb finite, 
but it still contains the IR divergences: 
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The real corrections arise from processes gg — > J/ijjgg, 
gg — > J/ipqq and gq{q) — > J/ipgq(q). The phase space 
integration of above processes will generate IR singular- 
ities, which are either soft or collinear and can be con- 
veniently isolated by slicing the phase space into differ- 
ent regions. We use the two-cutoff phase space slicing 
method [13], which introduces two small cutoffs to de- 
compose the phase space into three parts. And then the 
real cross section could be written as 



HC 



HC 



HC 
add ■ 



(6) 



The hard noncollinear part a HC is IR finite and is nu- 
merically computed using standard Monte-Carlo integra- 
tion techniques. a s from the soft regions contains soft 
singularities and is calculated analytically under soft ap- 
proximation. It is easy to find that soft singularities for 
a gluon emitted from the charm quark pair in the S- 
wave color singlet J/ip are canceled by each other. a HC 
from the hard collinear regions contains collinear singu- 
larities which is factorized and partly absorbed into the 
redefinition of the PDF (usually called mass factoriza- 
tion [23|). Here we adopt a scale dependent PDF using 
the MS convention given by [22j. After the redefinition 
of the PDF, an additional term cr^ d is separated out. 
Finally, all the IR singularities are canceled analytically 



for d"~ 



- r HC 



a v . And it is found that only one color 



factor d a i,c appears in both M B and M v with a, b and c 
being the color index of the three gluons in the process. 
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To obtain the transverse momentum distribution of 
J/tp, a transformation for integration variables (cLz^di — ► 
Jdptdy) is introduced. Thus we have 

= f JdxidyG g (xi,Hf)G g (x2,^f)^, (7) 

where y and pt is the rapidity and transverse momentum 
of J/tp in laboratory frame respectively. The polarization 
measurement a is defined as: 



da T /dp t - 2da L /dp t 
da T /dpt + 2da L /dp t 



(8) 



It represents the measurement of J/ip polarization as 
function of pt- To calculate a(pt), the polarization of J/tp 
must be kept in the calculation. The partonic differen- 
tial cross section with J / ip polarized could be expressed 
explicitly as: 



d(7A 
dt 



ae(X)-e*(X)+ £ oy Pi ■ e(A) Pj • e*(A), (9) 



where A = T U T 2 ,L. e(Ti), e(T 2 ), e(L) are the two 
transverse polarization vectors and the longitude one for 
J/tp, and the polarization of all the other particles are 
summed up in n-dimensions. It causes more difficult ten- 
sor reduction path than that with all the polarization 
being summed over in virtual correction calculation. It 
is founded that a and finite when the virtual cor- 

rection and real correction are summed up. 

For a^jc anc j fiS _j_ £HC _|_ frV ^ fo Q ca i cu i a ti on j s d one i n 

Eq. © as well as in the usual way in which the polariza- 
tion for all particles are summed up. The two results are 
used to check with each other numerically. In the third 
way, the polarization of gluon is also kept and used to 
check gauge invariance by replacing the gluon polariza- 
tion vector to the gluon's 4-momcntum in the final nu- 
merical calculation. To calculate a HC in real correction 
processes, the numerical amplitude calculation is used 
and the polarization of gluons are only summed for their 
physical freedom to avoid the involving of diagrams with 
external ghosts lines. 

In our numerical calculations, the CTEQ6L1 and 
CTEQ6M PDFs, and the corresponding fitted value for 
a s {M z ) = 0.130 and a s (M z ) = 0.118, are used for 
LO and NLO predictions respectively. For the charm 
quark mass and the wave function at the origin of J/tp, 
m c = 1.5 GeV and |i? s (0)| 2 = 0.810 GeV 3 are used. The 
two phase space cutoffs 8 S and 5 C are chosen as <5 S = 10~ 3 
and S c = S s /50 as default choice. To check the two cut- 
offs invariance for the final results, different values of 5 S 
and <5 C , down to 6 S = 10~ 5 , are used and the invariance 
is found within the error control (less than one percent). 
It is known that the perturbative expansion calculation 
is not applicable to the regions with small transverse mo- 
mentum and large rapidity of J/tp. Therefore, the result 
is restricted in the domain p t > 3 GeV and less 
than 3 or 0.6. 
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FIG. 2: Transverse momentum distribution of differential 
cross section with /v = fif = y (2m c ) 2 + at LHC (upper 
curves) and Tevatron (lower curves). Center mass energy are 

^Tcvatron = I- 98 TeV and \/«lhc = 14 TeV - NLO+ denotes 
result including contribution from gg — * J/tpcc at NLO. 



The dependence of the total cross section at the renor- 
malization scale \i r and factorization scale /if are ob- 
tained and it agrees with the Fig. 3 in ref [13 . In Fig.0 
the pt distribution of J / ip is shown. The pt distribution 
of J/tp polarization a is shown in Fig. [3] At LO, a is 
always positive and becomes closer to 1 as pt increases 
from 3 GeV to 50 GeV, which indicates that transverse 
polarization is always larger than longitude polarization 
and transverse polarization plays the major role in high 
Pt region. But there is extremely change when includ- 
ing NLO QCD corrections, a is always negative and 
becomes closer to -0.9 as pt increases from 3 GeV to 50 
GeV, which indicates that transverse polarization is al- 
ways smaller than longitude polarization and longitude 
polarization plays the major role in high p t region. Mean- 
while the J I polarization of process gg — > J/tpcc is near 
zero. By including contribution of this subprocess, the 
result shown in the figures as NLO + is closer to the ex- 
perimental result. 

In conclusion, we calculated the NLO QCD correction 
of J/tp hadronproduction at Tevatron and LHC. The 
method of dimensional regularization is taken to deal 
with the UV and IR singularities in the calculation, and 
the Coulomb singularity is isolated and absorbed into the 
cc bound state wave function. To deal with the soft and 
collinear singularities, the two-cutoff phase space slicing 
method is used in the calculation of real corrections. Af- 
ter adding all contribution together, a result which is 
UV, IR and Coulomb finite is obtained. Numerically, we 
obtain a K factor of total cross section (ratio of NLO 
to LO) of about 2 at fi r = fif = \J (2m c ) 2 +p\. The 
transverse momentum distribution of J/tp is presented 
and it shows that the NLO corrections would boost the 
differential cross section more and more as p t becomes 
larger and reaches about 2 or 3 order of magnitude at 
Pt = 50 GeV. It confirms the calculation in Ref. p7| . 
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FIG. 3: Transverse momentum distribution of polarization a with /x r = fj,f = (2m c ) 2 + at Tevatron(left) and LHC(right). 
The unlabeled dotted line denotes the polarization of gg — > J/ipcc and NLO + denotes result including that of gg —* J/ipcc. 



The real corrections for gg — > J/ipcc is also calculated 
and agree with those of Ref . [H[ . 

The J/ip polarization at NLO is studied for the first 
time and the results show that the J ftp polarization is 
extremely changed from more transversal polarization 
at LO into more longitudinal polarization at NLO. Al- 
though it gives more longitudinal polarization than the 
recent experimental result [l6| on the J/ip polarization 
at Tevatron. It sheds light on the solution to the large 
discrepancy of J/ip polarization between LO theoretical 



predication and the experimental measurement, and sug- 
gests that the next important step is to calculate the NLO 
correction for color octet state J/tp hadronproduction. 
To re-fix the color-octet matrix elements and to see what 
happens to the polarization of J/tp at NLO. 
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